Ionization of Rydberg atoms in a low frequency field: 
modelling by maps of transition to chaotic behavior 
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Abstract. We investigate a microwave ionization of highly excited atom in a 
low frequency field and show that such a process may be studied on the bases 
of map for the electron energy change during the period of the electron motion 
between two subsequent passages at aphelion. Simple approximate criterion 
results to the threshold field for transition to chaotic behavior very close to the 
numerical results. We show that transition from adiabatic to chaotic ioniza- 
tion mechanism takes place when the field frequency to the electron's Kepler 
frequency ration approximately equals 0.1. 



1. Introduction 

Highly excited atom in a monochromatic field is one of the simplest real strongly 
driven by an external driving field nonlinear system with stochastic behavior. That 
is why recently a large amount of effort, both theoretical and experimental, has been 
devoted to the studies of the Rydberg atoms in strong microwave fields (for review see 
[1-3] and references therein). The observation of excitation and ionization rates provides 
evidence for stochastic behavior of weakly bound electron: the ionzation of Rydberg 
atoms exhibits a threshold dependence on the field amplitude and, at least for the low 
relative frequency of the field, appears as a diffusion like process. The classical dynamics 
of the excited electron may be described by the map, called the "Kepler map", rather 
than by differential equations of motion [4,5], which greatly facilitates the numerical 
and analytical investigation of stochasticity and ionization process. The map, mostly 
represented for the number of absorbed photons, is widely used for the analysis of chaotic 
processes for relatively high frequencies of the field [1-5]. In the very low frequency region 
the ionization is direct: the electron escapes over the potential barrier [5-6]. However, the 
threshold field for the direct low frequency (adiabatic) ionization is considerably higher 
than that for the transition to chaotic behavior in the medium and high frequency field. 
On the other hand, the derivation of the mapping equations of motion [4,5] is based 
on the classical perturbation theory. Therefore, the question of the validity of these 
maps for description of the ionization process in the low frequency relatively strong field 
arises. In addition, the transition from adiabatic to chaotic ionization mechanism is of 
great interest, especially as the classical ionization theory is proper for the low frequency 
ionization. 

It is the purpose of this work to investigate the fitness of the classical maps [4,5] for 
the low-frequency ionization of Rydberg atoms and on this basis to analyse a transition 
from adiabatic to chaotic motion in the strongly driven by the slow external field system. 



2. Equations of motion and dynamics 




The direct way of coupling the electromagnetic field to the electron Hamiltonian is 
through the A • P interaction, where A is the vector potential of the field and P is the 
generalized momentum of the electron. Thus the Hamiltonian of the hydrogen atom in 
a linearly polarized field F cos(ut + where F and uj are the field strength amplitude 
and frequency, in atomic units has the form 
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The electron energy change due to the interaction with the external field follows from 
the Hamiltonian equations of motion [7] 

E — —r ■ Fcos(ujt + i?). (2) 

Measuring the time of the field action in the field periods one can introduce the scale 
transformation [8] when the scaled field strength and energy are F s = F/uo 4 / 3 , and 
E 8 = E/u 2 / 3 . However, it is convenient to introduce the positive scaled energy e = —2E S 
and relative field strength F = Fn 4 = F s /eq, with n being the initial principal quantum 
number. 

Note, that eq. (2) is exact if f is obtained from the equations of motion including the 
influence of the electromagnetic field. Using the parametric equations of motion we can 
calculate the change of the electron's energy in the classical perturbation theory approx- 
imation. We restrict our subsequent consideration to the one-dimentional model, which 
corresponds to the states with low orbital quantum numbers / <C n. The integration of 
eq. (2) for the motion between two subsequent passages at the aphelion (where x — 
and there is no energy change) results to the map (see [4,5] for details) 
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where s = e~ 3 / 2 = 2E) 3 / 2 is the relative frequency of the field, J' s (z) is the derivative 
of the Anger function and function G may be obtained from the requirement of area- 
preserving [4]. In the low frequency s C 1 limit J' s (s) ~ s/2 and we have the map 




- 2nF Q (el/e%) sintf, 
■&j + 2-Kje% + 5ttF (e 2 /e%) cos^ 



(4) 



The results of the numerical analysis of maps (3) and (4) in the low frequency, s < 1, 
limit are presented in Fig. 1 and 2. We see that the threshold ionization field approaches 
the static field ionization threshold F st 0.13 when s — > (e — > oo). Thus, the maps 
(3) and (4) are valid in the low frequency limit when the strength of the driving electric 
field is of the order of the Coloumb field. 

For the low frequencies, 2irs = 2ti je 3 ^ 2 1, the change of the angle d after one step 
of iteration is small and we can transform the system of equations (4) to the differential 



Fig. 2. The relative threshold field strengths for the ionization outset from the numerical 
analysis of the maps (3) and (4) and according to the approximate criterion (8)-(9) 



equation 

rf(cos^) e 5cos$ . . 

de = ejF~o + 2e ' ^ 
The analytical solution of eq. (5) with the initial conditions e = e when •§ = $ is 



costi = z 5 cos$ - (2/F )z 4 (l- z), z = Je/7 (6) 

For the relatively low values of F there is a motion in all the interval 4- 2tt of the 
angle -d. However for F = 2z 4 /5 the increase of ■& at d ~ 7r changes to the decrease and 
results to the fast decrease of e and ionization process. The minimal value of F for such 
a motion corresponds to $o = and may be defined as a maximal value of F resulting 
to the motion in the interval 4- 2it, i. e., the maximum of the expression 

F = 2z A (l-z)/(l + z 5 ). (7) 

Such a maximum is at z = z , where z is a solution of equation z 5 + 5z — 4 = 0, i. 
e., z = 0.75193 and results to the F ° = 0.1279 which is only 1% lower the adiabatic 
ionization value F S * = 2 10 /(3tt) 4 ~ 0.1298. 

The approximate criterion for transition to chaotic behavior is [5] 

K = max IS&j+JS&j - 1 1 > 1 (8) 

and according to eqs. (4) yields to the expression for the threshold field strength 

F C = £ 5 /6tt 2 ^ = zl°/6n 2 sl (9) 



where z c is the solution of eq. (7) with F = F C . Note that zf ~ 1 - 10F C + 30(F C ) 2 =f . . . 
if Fq < 0.1. For 0.09 < s < 0.5 expression (9) gives the ionization threshold field very 
close to the numerical results (see Fig. 2). 

3. Conclusion 

From the above analysis we can conclude that the map at aphelion (3) is suitable for 
the investigation of transition to chaotic behavior and ionization of Rydberg atoms also 
in the low frequency field, even for adiabatic ionization, when the strength of the external 
field is comparable with the Coloumb field. For such a purpose there is no need to use 
the map for the two halves of intrinsic period [9]. Moreover, the approximate criterion 
(8) for transition to chaotic behavior yields to the threshold field strength very close to 
the numerical results if we take into account the increase of the electron's energy by the 
influence of the electromagnetic field. The transition from adiabatic to chaotic ionization 
occurs at the relative field frequency s ~ 0.1. 
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Fig. 1. Trajectories of the map (3) for the different initial conditions Eo, 0o and relative field 
strength F . The pictures in the left-hand side correspond to the regular quasiperiodic motion, 
while in the right-hand side represent the ionization process for a little bit stronger field. At 
£o=4.3, i. e. at 50=0.11 the transition from the adiabatic to the chaotic ionization mechanism 
takes place. 



